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Abstract: We compute the virtual O (a^) fermionic contributions to the charm quark 
mass dependent matrix elements of the B Xg'f decay. In the case of a massless quark 
loop insertion into the gluon propagator, our result obtained as an expansion in z = m^/m^ 
and an exact expression in terms of multi-fold MB integrals, confirms the findings of Bieri, 



Greub and Steinhauser |2C]. We observe, however, large deviations in the case of a b-quark 



loop insertion. The charm quark loop shows smaller, but still noticeable differences. 
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1. Introduction 

One of the most interesting rare B-meson decays is the inclusive B — > Xg'j mode which 
provides precise and clean short-distance information on AB = 1 flavor changing neu- 
tral currents (FCNCs). This process occurs in the SM only at the loop level, through 
the exchange of W-bosons and up-type quarks, making it highly sensitive to non-standard 
effects which are not suppressed by additional factors a relative to the SM contributions. 
Combined with the low sensitivity of B ^ Xgj to non-perturbative effects, this makes it 
possible to observe new physics contributions indirectly, or to set limits on the relevant 
masses and coupling parameters. In fact the decay width F — >■ Xsj) is well approx- 
imated by the partonic decay rate T(b Xf"^*°"'7) which can be analyzed within the 
framework of renormalization-group-improved perturbation theory . 

In view of the above, it is obvious that both accurate measurements and precise theo- 
retical SM calculations, with a good control of both perturbative and non-perturbative cor- 
rections, have to be provided. As far as the experimental side is concerned, measurements 
are performed by several B physics experiments operating within various experimental set- 
tings: CLEO § (Cornell), BaBaR § (SLAC), Belle Q (KEK) and ALEPH |] (CERN) . 
Combining the measurements of the first three (the last one has very large error bars) for 
the branching ratio B{B Xg'j) leads to a world average with a cut > 1.6 GeV in the 
-B-meson rest frame which reads 

B{B ^ X,7)i;>,aGcv = (3.55 ± 0.24 ± 0.03) x 10-^ (1.1) 

where the first error is given by the statistic and systematic uncertainty, the second one is 
due to the theory input on the shape function, and the third one is caused by the 6 — > (^7 
contamination . This average is in good agreement with the recent theoretical estimate 
including known next-to-next-to-leading-order (NNLO) effects [^j 

B{B ^ X,7)r>i.6Gcv = (3.15 ± 0.23) x 10-^ (1.2) 
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where the error consists of four types of uncertainties added in quadrature: non-perturbative 
(5%), parametric (3%), higher-order (3%) and mc-interpolation ambiguity (3%). 

The total experimental error of about 7% in Eq. is of the same size as the 



expected 0{al) corrections to the perturbative transition b Xf^^^°^j [^], which calls for 
completing the SM calculations with this accuracy level. 

QCD corrections to the partonic decay rate F (6 — > 57) contain large logarithms of the 
form a" (m^) In™ (nih/M), where M = mj or M = myy and m < n (with n = 0, 1, 2, . . . ) . 
In order to get a reasonable prediction for the decay rate with next-to-next-leading-log 
(NNLL) precision, it is necessary to resum logarithms with (m = n ,n — 1, n — 2) with 
the help of renormalization-group techniques. A convenient framework is an effective low- 
energy theory obtained from the SM by decoupling the heavy electroweak bosons and the 
top quark. The resulting effective Lagrangian, given in the next section, is a product of 
the Wilson coefficients Cj(/i) which play the role of coupling constants, and local flavor- 
changing operators Qi{n) ■ 

As far as the next-to-leading order corrections are concerned, the program has been 
completed already a few years ago, thanks to the joint effort of many groups (see e.g. |9| 
and references therein). Moreover, each of the ingredients has been cross-checked by more 
than one group^. The next-to-next-to-leading order calculation, which involves hundreds 
of three-loop on-shell vertex-diagrams and thousands of four- loop tadpole-diagrams, is a 
very complicated task and is currently under way . A consistent calculation of 6 — > 57 at 
this order requires three steps: 

• Matching: Evaluation of Cj(^o) ^-t the renormalization scale ^0 ~ Mw,mt by re- 
quiring equality of the SM and effective theory Green's functions at the leading order 
in (external momenta) / {M\y , mt) to O (a^) . All the relevant Wilson coefficients 



have already been calculated |11, 12 1 to this precision, by matching the four-quark 
operators Qi, ■ ■ ■ ,Qq and the dipole operators Qj and Qs at the two- and three-loop 
level respectively. 

Mixing: Calculation of the operator mixing under renormalization, by deriving the 
effective theory Renormalization Group Equations (RGE) and evolving Cj(^) from ^0 
down to the low-energy scale fib ~ using the anomalous-dimension matrix (ADM) 
to O (af) . Here, the three-loop renormalization in the {Qi, . . . , Qe} and {Qj, Qs} 
sectors was found in |]l3|, and results for the four-loop mixing of Qi, . . . , Qe into 



Qy and Qs were recently provided in |15], thus completing the anomalous-dimension 
matrix. 



• Matrix elements: Evaluation of the on-shell 6 ^ 57 amplitudes at fib ~ mb to 
O (ct^). This task is not complete yet, although a number of contributions is known. 
The two- loop matrix element of the photonic dipole operator Qj, together with the 

^In fact the complete calculation of Bremsstrahlung corrections has not been checked, but the effects 
are extremely small. 



- 2 - 



corresponding bremsstrahlung, was found in 17 1, confirmed in [18| and subse- 
quently extended to include the full charm quark mass dependence in [|^]. In pO| , 
the O (a^nj) contributions were found to the two- loop matrix elements of Qj and 
Qs, as well as to the three- loop matrix elements of Qi and Q2, using an expansion 
in the quark mass ratio ml/m^. Diagrammatically, these parts are generated by in- 
serting a one-loop quark bubble into the gluon propagator of the two-loop Feynman 
diagrams. Naive non-abelianization (NNA) is then used to get an estimate of the 
complete corrections of O (a^) by replacing Uf with — |/3o. Moreover, the contri- 
butions of the dominant operators at O (a^/^o) to the photon energy spectrum have 



been computed in |21]. Finally, in [^], the full matrix elements of Qi and Q2 have 
been computed in the large rric limit, rUc ^ mi,/ 2. Subsequently, an interpolation in 
the charm quark mass has been done down to the physical region, under the assump- 
tion that the /3o-part is a good approximation at nic = . This is the source of the 
interpolation uncertainty that has been mentioned below Eq. (|L 



To date, no independent check has been provided for the O (a^nj) results for the ttt-c- 
dependent matrix elements of Qi and Q2, despite the fact that they constitute a major input 
both for the NNA and for the interpolation of the non-NNA terms between rUc ^ and 



rUc < 'mb/2, and are thus crucial for the accuracy of Eq. (1.2). Besides, for the numerical 
evaluation, it was assumed that nj = 5 massless fermions are present in the quark loop. 
Since the charm and bottom quark masses are not negligible, it is important to check their 
numerical relevance. 

In this paper, we present a calculation of the virtual O (a^n/) contribution to the 
matrix elements of the operators Qi and Q2- As far as contributions from diagrams with 
a massless quark loop insertion into gluons are concerned, our result was obtained using 
two different methods: an expansion in the mass ratio m^/m^, which confirms the findings 



of |2C], as well as an exact evaluation in terms of multi-fold Mellin Barnes (MB) numerical 
integrals. The complexity of diagrams with a massive loop insertion (charm and bottom) 
required a mixed approach, where both MB techniques and differential equations were used 
numerically. 

This paper is organized as follows . In the next section, we introduce the relevant 
effective Lagrangian and present the methods used. Subsequently, our results for the matrix 
element of the operator Q2 {Qi is related to Q2 by a color factor at this level of perturbation 
theory) are given in the form of fitting formulae, which satisfactorily approximate both the 
massless and the massive cases over the full experimentally interesting range of values of 
the charm and bottom quark mass ratio. Finally, we give our conclusions. An appendix 
contains the behavior of the contributions in the limit nic ^ mb- 

2. Calculation 

As mentioned in the introduction, we work within an effective low-energy theory with five 
active quarks, obtained from the SM by integrating out the heavy degrees of freedom with 
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a mass scale M > M^. The Lagrangian of such a theory reads 

Aff = -CQCDxQEDltt, d, s, c, b) + —^VtlVtb 2^ Ciiij) Qi{ii). (2.1) 

Here the first term is the usual QCD-QED Lagrangian for the light SM quarks. In the 
second term, Vij denotes the elements of the Cabbibo-Kobayashi-Maskawa matrix, Gf 
is the Fermi coupling constant and Ci{^) are the Wilson coefficients of the corresponding 
operators Qi evaluated at the scale ^. Adopting the operator definitions of ||2^ , the physical 
operators that are relevant for our calculation read 



Q2 = (sL7/iCL)(cL7^6L), 

Q4 = {sli^t%l)Y1 



Qt = -^mbiti)isLaf''bR)F^,, (2.2) 



where Qi and Qi are needed for the renormalization, and rn^(/x) is the b-quark MS mass. 
The sum over q runs over all light quark fields, and e and gs are the electromagnetic and 
strong coupling constants respectively. -P^jy is the electromagnetic field strength tensor, 
and (a = 1 ... 8) denote the SU{3) color generators. 

The choice of the basis of the local four-quark operators is made such that no problems 
connected with the treatment of 75 in d = 4 — 2e arise [23|. This allows for consistent use of 
a fully anticommuting 75 in dimensional regularization throughout the whole calculation. 
The only evanescent operator needed for renormalization would in principle be Qn defined 
as 

Qn = {sl1^.1.1pT''cl){cl1^YYT%l) - 16 Qi , (2.3) 

but it turns out that it does not contribute to the fermionic part. 

Since one can show that at this order, the contributions to the matrix elements of the 
operator Qi are related to those of Q2 by a simple substitution of color factors 

{si\Qi\b) = -^^{sj\Q2\b) (2.4) 

we concentrate in the following on the derivation of the result for Q2. Upon reducing the 
set of scalar integrals involved in the calculation of diagrams with a Q2 insertion, of which 
a sample is shown in FIG. |l], we are left with a set of 18 master integrals in the case with 
a massless quark loop insertion into the gluon propagator, 47 master integrals in the case 
with a massive b-quark loop insertion and 38 for the c-quark case. These master integrals 
have been calculated by combining two different approaches. 

In the first approach, we have used the Mellin-Barnes technique |2^, which relies 
on the identity 

1 _ r^+*°° r(A + z)r(-z) dz 

WTW'Jp-ioo nx) 2^' 



(2.5) 
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Figure 1: Typical 3-loop vertex graphs calculated in this paper. The dotted bubble refers to a massless 
quark. 



where —Re A < /3 < 0. With this relation a sum of terms raised to some power is replaced 
with a product of factors. We have derived our MB representations for the master integrals 



using the automated package |2G|, and used the MB package of ref. [27| for their analytic 
continuation. In the case where the inserted quark loop into the gluon propagator is 
massless, we have used two methods. In the first, we have performed an expansion in the 
quark mass ratio z = m'^/m'^ by closing contours. The coefficients of this expansion were 
given by (at most) one-dimensional MB integrals which have subsequently been expressed 
as a sum over residues and resummed with the help of XSummer pS]. Our result for the 



matrix element of Q2 using this method is in full agreement with |20|. It is also consistent 
with the second method, where the exact z-dependence was retained and a numerical 
integration of the MB representations was performed using the MB package. 

Due to poor convergence, it turned out not to be possible to compute all the master 
integrals that occur in diagrams with a massive quark loop insertion with the help of MB 
representations. For these particular cases, we have used a second approach based on the 
method of differential equations. Using the fact that the master integrals Vi{z,e) (after 
rescaling by a trivial factor) are functions of e and the mass ratio z or its inverse y = z~^, 
respectively, a system of differential equations has been generated, 

^Vi{y,e) = Ai,j{y,e)Vj{y,e), (2.6) 

where the right-hand side was again expressed through master integrals with the help of 
relations obtained from the reduction. The block-triangular matrix Aij{y,e) is composed 
of rational functions of e and y. 

The solution of this system for arbitrary values of y proceeded in two steps. First, an 
expansion in £ and y for e, y — > was performed with the ansatz 

y,(?/,e) = ^Ci„^fce"y'"log'=y, n, m = -3, -2, . . . ; = 0, . . . , 3 + n G(n) (2.7) 

nmk 



and the coefficients were calculated recursively up to high powers of y |29|. In this limit, 
m'^ ^ and the initial conditions correspond to vertex diagrams that can be derived in 
an automated way from diagrammatic large-mass expansions. With this series at hand we 



- 5 - 



0.04 



0.08 



0.12 



0.16 



0.2 



Figure 2: Plot of Re(s7|(52|^)ri^/ '° for one massless flavor and /i;, — nib- The dashed vertical line 
corresponds to the central value ml/ml — 0.068. 



were able to compute the master integrals for y <^ 1 with high precision. In the second 
step a numerical integration of the e-expanded system into the physical region, y > 1, was 
carried out using the high precision values at a starting point y <C 1. For this task, we 
have used the Fortran package ODEPACK |30|, which allows for numerical solutions of 
huge systems of differential equations. To avoid points of numerical instability on the real 
axis, the integration path has been shifted into the complex plane. 

Besides the bare three- loop matrix elements, we also needed matrix elements of two 
types of counterterms: those coming from the renormalization of the strong coupling con- 
stant as (known from the QCD beta function pT| , |3^), and those due to the mixing of 
Q2 into other operators. For the corresponding details, we refer the reader to ||2^. We 
just mention here that the renormalization constants used for the renormalization of the 
diagrams with a massless quark loop insertion into the gluon propagator are the same as 
those used for the renormalization of the diagrams with a massive quark loop insertion, 
because all the needed renormalization constants were derived in the MS-scheme and are 
therefore mass independent. 



3. Results 



We now present the results for the amplitudes of the three different quark loop insertions 
with the following normalization 

{si\Q2\b)o{a2nf) = (^) -^^b Uf {s-f\Q2\b)^n}'^^ UsR^^Ub (3.1) 

where nif, denotes the b-quark pole mass, e and q are the photon polarization and momen- 
tum, R = (l-|-75)/2 is the right handed projection operator, and nj is the number of active 
flavors of a given mass. The superscript (2) counts the powers of Us and M = (0, mhOicmc) 
denotes the mass of the quark running in the loop inserted into the gluon propagator. 

It is important that there is no need for a b-quark mass renormalization in our cal- 
culation, so we use the different schemes as guided by the complete calculation, and thus 
turn to the IS mass [^] for our final study. 



Input parameter 


experimental value 


nicimc) 


(4.68 ± 0.03) GeV || 

(1.224 ± 0.017 ± 0.054) GeV |35] 



Table 1: Experimental inputs relevant for the present calculation 



All our subsequent results are given in the form of fitting formulae that cover the whole 
interesting range of variation of z. With the current input given in Table ||, and allowing 
for a 3 sigma variation of mc{mf.) and m^^ , this corresponds to z G [0.017,0.155]. In the 
plots of Fig. ^ and Fig. |3|, we also show the central value, which is currently z = 0.068. 

In the case of a massless quark loop insertion we agree with the previous calculation 
pO|] . Nevertheless we present here a fitting formula for easier comparison with our new 
results 

Re(s7|Q2|&)^5'° = 9.080 - 0.7624 z - 5.069 + 12.61 z Inz 

+ (-9.679 + 5.157 z + 1.726 z^ - 16.18 z In z) ln(m6/^) 

+ ^ \n\m,/^,) (3.2) 

In Fig. ^ we show the data points used to obtain the fit function together with the result 
Eq. |3.2|. Motivated by the presence of logarithms in the small z expansion, we have included 



terms proportional to z\\iz above, and in equations (|3.3| ) and (3.4) below, which improved 
the quality of the fit significantly. As it stands, our fit function reproduces the exact values 
with a relative precision of at least 10~^. 

Our result for the contribution of the diagrams with a massive b-quark loop insertion 
is given by the fitting formula of the same relative precision as above 

Re(s7|Q2|^>^5'"'' = -1.836 + 2.608 2 + 0.8271^2 -2.441 z In z 

+ (-9.595 + 5.157Z + 1.726^2 -16.18Z \nz)\n{mb/ fi) 

+ ^ ^^^\m,/^,) (3.3) 

which is plotted in Fig. ^(a). It is interesting to note that the massless approximation 
overestimates the massive result by a factor of 6 and, moreover, has the opposite sign. This 
points to an expected decoupling like effect. Remark, however, that the scale dependence 
is just as strong as in the massless case. 

When the massless approximation is confronted with the result for the charm quark 
loop insertion, one still observes noticeable differences as shown in Fig. ^(b). The obtained 
fitting formula of relative precision of 10~^ is 

Re(s7|Q2|i>i^^^'"'" = 9.099 + 13.20 z - 19.68 + 25.71 z Inz 

+ (-9.679 + 13.62 z - 13.94 z^ - 12.98 z In z) ln(mb//i) 

+ ^ \n\m,/^,) (3.4) 
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Figure 3: Plots of Re(s7|Q2|fe)rf/ as function of m^/m^ with M 
fib — TTii,. For comparison, we also show the M = case. 
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When z approaches both Eq. 3^ and Eq. 3^ tend approxhnately to the same value, as 
it should be. The residual difference is due to the quality of the fit outside the validity 
range, and the steepness of the curves in Fig. ^ 

Let us finally comment on the renormalization group logarithms in Eqs. |3.2| to 3.4. It 



has already been shown in the massless case in [22| that the single logs can be expressed 



through the a{z) and b{z) functions known from the NLO calculation |36]. In fact the 
exact expression for the coefficient of log(m5//i) in this case is 



2Qn 

^ 5R(a(z)+6(z))- — 



Using the RGE one can show that the coefficient of the same log in Eq. 3.3 is 



and similarly for Eq. 3.4 



U^{a{z) + b{z)+b{l))-^ 



-i^^(a{z)+2b{z))-—y 



(3.5) 



(3.6) 



(3.7) 



4. Conclusions 

In this paper, we have presented a calculation of the 0{a1nf) contribution to the charm 
quark mass dependent matrix elements of the B — > Xg'j decay. In |20|, a similar calcu- 
lation was done as a mass expansion in the ratio z = m^/m^ assuming that the b and 
c quarks inside the fermionic loop inserted into the gluon propagator are massless. The 
results have then been used to estimate the complete corrections of O(a^) to the matrix 
elements {sj\Qi^2\b) by replacing nj with — |/3o) i-e. following the naive non-abelianization 
procedure. They form an important part of the NNLO contributions used in the recent 
theoretical estimate of the branching ratio B{B Xg'j) [^. Our goals were firstly to pro- 
vide an independent check of the O(a^nj) corrections for massless quarks, and secondly 



■^We thank M. Misiak for pointing this out. 
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to check the vahdity of the massless approximation by keeping the full mass dependence. 
These goals have been reached, and our result in the massless approximation, which we 
have obtained as a mass expansion in z as well as an exact evaluation in terms of multi-fold 
MB integrals, confirms the findings of |^^. However, although the massless approximation 
reproduces the contribution of the diagrams with a massive c-quark loop insertion reason- 
ably well, it is not justified for the diagrams with a massive b-quark. Neglecting the mass 
of this particle leads to an estimate of its contribution that is 6 times larger than its true 
value and moreover has the opposite sign. Of course, since this effect comes from a single 
quark family its impact on the branching ratio is rather mild. It turns out that, when 



compared with the interpolation from [ 22 1 , one finds an enhancement between one and two 



percent, depending on the renormalization scale [i. 
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A. Contributions in the limit ^ 

During our calculation we also derived the large limit of all the contributions. In the 



massless case, we agree with |22]. For the two massive cases, we have (with /i = mb) 



'\Sl\Q2\hfn}'^'' = 4.25648 0.503085 Inz 0.888889 In^ z 



i (-0.725053 - 1.80916 \uz + 0.0938272 In^ z) 

+ 4t (-1.39486 -0.968501 In z -0.147443 In^z) -hO ( 4 
z^ V 



(A.l) 



and 



(s7|Q2|^>E^'"''^ = 1.67932 ^ 0.526749 Inz 0.823045 In^ z 



^(0.20839 -^0.11775 lnz + 0.128395 \^ z) 

4 (-0.0360638 - 0.0470166 Inz + 0.0324515 In^ z) O f 4r 

z^ V z-^ 
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